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INTRODUCTION 


This  report  deals  with  the  mathematical  aspects  of  a problem  of 
importance  in  the  failure  of  structures,  namely,  the  sudden  spreading 
of  a crack  in  a stressed  medium.  A previous  study  on  the  dynamic  prob- 
lem of  torsion  of  a circular  crack  was  made  by  Sih  and  Embley  [1]. 

Their  work,  however,  is  restricted  to  the  case  of  a crack  which  opens 
instantaneously  to  a fixed  radius  and  thereafter  remains  stationary. 

Most  work  on  crack  propagation  under  other  loading  conditions  is  based 
on  the  mathematically  expedient  supposition  that  propagation  occurs  at 
constant  speed.  It  is  more  realistic  to  consider  the  speed  of  propaga- 
tion as  a function  of  time  based  on  some  physical  hypothesis.  Recently, 
Kostrov  [2]  has  proposed  an  analytical  technique  for  treating  a vari- 
able speed  crack  and  has  obtained  results  for  the  case  of  anti-plane 
loading  conditions.  In  this  report  Kostrov's  technique  is  applied  to 
the  case  of  to’rsional  loading  conditions. 

Kostrov's  method  has  the  distinct  advantage  over  other  methods, 
notably  the  Weiner-Hopf  method  and  the  dual  integral  equation  method, 
in  that  mixed  problems  with  time  dependent  boundaries  can  be  treated. 
Moreover,  the  near-field  solution  can  be  deduced  rather  easily,  at 
least  for  torsion  and  antirplane  strain  conditions.  An  inherent  draw- 
back of  Kostrov's  method  is  that  for  a finite  length  crack  the 

^G.  C.  SIH  and  G.  T.  EMBLEY  (1972)  Journal  of  Applied  Mechanics, 
395-400.  Sudden  Twisting  of  a Penny-Shaped  Crack. 

^B.  V.  KOSTROV  (1966)  Journal  of  Applied  Mathematics  and  Mechanics, 
30(6),  1241-1248.  Unsteady  Propagation  of  Longitudinal  Shear 
Cracks. 
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analytical  results  are  limited  to  the  time  interval  for  which  waves  ' 
enamating  from  opposite  points  on  the  crack  border  have  not  yet  come 
into  contact;  it  appears  that  numerical  integration  is  required  for  the 
time  in  which  wave  interaction  occurs.  On  the  basis  of  the  results  of 
the  study  by  Sih  and  Embley  [1],  however,  the  dynamic  stress  intensity 
factor  reaches  its  maximum  value  before  wave  interaction  occurs.  Thus 
it  is  reasonable  to  suppose  that  the  results  obtained  in  this  report  on 
the  dynamic  stress  intensity  factor  represent  the'maximum  value  that 
can  be  expected. 

PROBLEM  FORMULATION  AND  SOLUTION 

Consider  an  infinite  elastic  solid  with  shear  modulus  y and  den- 
sity p referred  to  cylindrical  coordinates  r,0,z.  At  time  t=o  a cir- 
cular crack  of  radius  a occupies  the  plane  z=o  ; for  time  t>o  the 
crack  radius  1.s  r*(t)  = a + R(t),  where  R(t)  is  a monotonic  increasing 
function,  with  R(o)=o.  It  is  assumed  that  ^<c,  where  c=(y/p)^^^ 
is  the  torsional  wave  speed  of  the  solid.  Suppose  that  the  solid  is 
loaded  by  tractions  of  equal  magnitude  applied  to  the  adjacent  edges 
of  the  crack  surface,  and  that  these  tractions  are  independent  of  0. 
Then  the  displacement  v(r,z,t)  will  be  anti -symmetric  in  z,  so  that  it 
is  sufficient  to  consider  only  the  region  z>o.  For  this  problem  the 
stress  components  are 

’g.  C.  sih  and  G.  T.  EMBLEY  (1972)  Journal  of  Applied  Mechanics. 

395-400.  Sudden  Twisting  of  a Penny-Shaped  Crack. 
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(2) 


and  the  displacement  satisfies  the  equation 

\rr  " ; \r  - 72  " \zz  = \tt  * 

wherein  for  convenience  ct  has  been  replaced  by  t.  The  boundary  con- 
ditions on  z=o  may  be  written  as 

(3) 

V = 0 , r>r*(t)  , (4) 

where  f(r,t)  is  a continuous  function.  Since  by  superposition  any 
initial  displacements  may  be  separated,  then  without  loss  in  generality 
let 


Ogz  = f(r,t)  , r<r*(t) 


V = 0,  V . = 0 

9 ^ 


t = 0 . 


(5) 


The  problem  as  formulated  by  equations  (l)-(5)  constitutes  a mixed 
boundary  value  problem.  In  order  to  apply  Kostrov's  method  it  is 
first  necessary  to  obtain  a relationship  between  the  displacement  and 
the  stress  on  the  plane  z=o.  To  this  end  let 


v(C,z,p)  = ^e’P^dt  /“  rJ^(Cr)  v(r,z,t)  dr  , (6) 

E(^,z,p)  = J^e'P^dt  /“  rJ^(^r)  a(r,z,t)  dr  , (7) 


denote  Hankel -Laplace  transforms;  in  (7)  a denotes  a.  . Applying 

oZ 

(6)  to  equation  (2)  and  using  (5)  gives 

O 

v^2z  ■ “ ^ ° 

where  c^=  + p^.  The  appropriate  solution  of  (8)  is 

v(^,z,p)  = v(^,p)  , (9) 
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provided  that  Re{a}>o.  Likewise,  applying  (7)  to  (1)  and  using 
(9)  gives 

E(C.z,p)  = E(C,p)e”“^  , (10) 

and 

= - yav(5,p)  . (11 ) 

Formula  (11)  gives  the  relationship  between  the  transformed  stress  and 
displacement  components  on  z=o.  Now  substituting  (11)  into  (9)  and 
applying  the  Hankel -Laplace  Inversion  formula  gives 

pv(r,z,t)  = - ^ ASICS'") 

where  B denotes  the  well-known  Bromwich  contour.  Using  the  result 


1 f pt  . 

= 


a 


o<t<z 


Jq(5u)  , z<t<®  , 


(13) 


where  u = (t^-z^)^^^  , together  with  the  convolution  theorem  gives 


yv(r,z,t)  = - dt'  /“  r'a(r'.t')  dr'  (r^)J^  (r ' ^)J^(s^)d5, 

(14) 

2 2 1/2 

wherein  s = [(t-t‘)  -z  ] ' . The  inner  integral  may  be  expressed 

as  [3] 


/”cJi(r?)Jl(r'?)Jo(sC)  d5  = 


0 , o<r  < s-r  , 


Trrr 


; cot(})  , |s-r|<r'<s+r  , 


0 , r+s<r'  , 


(15) 


‘I.  S.  GRADSHTEYN  and  I.  M.  RYZHIK  (1965)  Tables  of  Integrals, 
Series,  and  Products,  Academic  Press. 
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where 


cot({)  = (r^+r'^-s^)  {[s^-(r-r‘ )^]  [(r+r‘ . (16) 

Hence  (12)  becomes 


yv(r,z,t)  = - o(r',t')  cot4>  dr'dt'  . 


(17) 


Region  R is  depicted  in  figures  (1)  and  (2)  for  the  situations  in  which 
o<t<(r^+z^)^ and  (r^+z^)^^^  <t,  respectively.  In  figure  (2)  region 
R is  bounded  by  the  line  t'=o  and  the  two  hyperbolas  s=r+r'  and 
s=|r-r' I . The  dashed  lines  represent  the  limiting  boundary  of  region 
R as  z-K),  and  in  the  sequel  will  be  referred  to  as  region  R^. 

Formula  (17)  may  be  used  in  conjunction  with  Kostrov's  method  to 
determine  the  stress  ahead  of  the  crack  perimeter.  Thus  putting  z=o 
in  (17)  and  invoking  conditions  (3)  and  (4)  gives 


/ / o(r',t')  cotcf)  dr'dt'  = - / / f(r',t')  cot4>  dr'dt',  r>r*(t) 


R 


o2 


R, 


ol 


(18) 


where  4>  = » and  where,  as  shown  in  figure  (3)  regions  R^^  and 

R^2  denote  portions  of  region  R^  to  the  left  and  right,  respectively, 
of  the  crack  perimeter  r'  = r*(t').  It  is  clear  that  o(r',t')=o  in 
region  R^^  • 

Equation  (18)  is  now. recast  in  terms  of  the  coordinates 
= r'-t'  and  n'  = r'+t'.  Thus  using  the  notation  a(r',t')  = 
T(?',n')  . f(r',t')  = g(?',n')  gives 


5 
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Figure  2.  Solution  Region:  Drawn  for  the  Case  t>R,  r>r*(t) 


V 

+ 

i. 
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Figure  3.  Diagram  for  Crack  Plane  Stress  Solution 


x(^',n') 

n*(^') 


N(g' ,g,n‘ »n) 

(n-n')’^^(cn')’'^ 


dn' 


* F(5'.E.r) 


where 


F(r,5,n)  = g(5'.n’)  — dn'  , 

5'  (n-n')’/2(5t^.)l/2 


(19) 


(20) 


N(5',^,n',n)  = 2 (?'n'nn)  + (n'-^‘)  (n-?)  , (21) 

and  n'=n*(?‘)  denotes  the  crack  perimeter.  Note  that  since  r'=r*(t') 
then  t'=r*"^(r')  = t*(r'),  so  that  n*(^')  is  given  implicitly  by  the 
formula 

n*  = ^'  + 2 t*  [ ^(n*+5‘ ) ] . (22) 

Now  equation  (19)  is  satisfied  if  the  expression  within  braces  vanishes. 
Moreover  substituting 

t(5'.h')  = (5+n')’^S(5'.5.n')  . (23) 

leads  to  an  Abel  type  integral  equation  for  the  stress  ahead  of  the 
crack  perimeter  of  the  form 


;;(^,)MC'.^.n') 


(n-n')''^  (2C-5'tn')  ♦ 


(g+n')(g'-^n') 

(n-r')'/2 


dn'  = 


-F(5'.5,n)  . . (24) 

The  solution  of  (24)  may  be  obtained  in  a straightforward  way  by  use 
of  Laplace  transforms;  the  result  is  written  as 


r [(5'+n)^(5+n)’'^^  t(5’,5.ri)] 


F(5'.5 

(n-n' ) 


dn'  . 


i 9+^'  d^  ,n 

^ (nH)^/^  dn^  n*(^') 


(25) 


9 


To  get  the  solution  at  the  point  (^,ri)  let  and  substitute  frorn 
(23).  This  gives 


y [(n+5)^  T(c.n)]  = - - 
dri  7T 


n*(5) 


F(C>g.n' ) 
(n-n' 


dn'.  (26) 


Now  since  F(5,5,n')  = (^+n‘)  E(^,n'),  where 

1 /2 

E(C.n')  = J*^*^^*  9(t.v)  dv  , (27) 

then  (26)  may  be  written  as 


j-  [(nn)2T(5.n)]  = - 1 f [(n+5)^'^  /\,  ,^'•^’’''1,,  dn'].(28) 

dn  TT  dn  dn  n*(E)(n-n')^/2 


Finally,  substituting  from  (27)  and  performing  one  integration  yields 
the  result 
x(E.n)  = 


1 


n-n' 


ii(rin)'/2[n-r*(5)]'/2  ? 

dn'  . n>n*(E)  • 


g(5.n') 


In  terms  of  the  original  variables  (29)  becomes 

1 


a(r,t)  = - 


^rV2(-^_r*(^)]l/2  |r-t| 


»t-r+r ' ) 


(29) 


r'^^^[r*(T)-r']^^^ 

r-r' 


dr' 


r>r*(t)  . 


(30) 


In  formula  (30)  the  quantity  t,  (not  to  be  confused  with  the  stress 
T(E.n))>  is  given  implicitly  by  the  relationship  t-T  = r-r*(T),  as  may 
be  seen  from  figure  (3).  The  quantity  t represents  the  retarded  time 
associated  with  the  moving  crack  border,  that  is,  a point  r experiences 
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at  time  t the  disturbance,  in  the  form  of  a wave,  caused  by  the  moving 
crack  border  when  the  border  location  was  at  r*(x). 

The  result  (30)  is  valid  for  the  time  interval  r-a<t<r+a.  This 
interval  corresponds  to  the  time  when  the  disturbance  from  the  moving 
crack  perimeter  point  (r*(t),  6+Tr)  has  not  yet  reached  the  crack  posi- 
tion (r*(t),6).  To  determine  a(r,t)  for  larger  values  of  time  it 
appears  necessary  to  replace  f(r',t')  by  -a(r',t')  on  those  portions  of 
the  range  of  integration  on  which  the  stress  is  unknown.  This  proce- 
dure corresponds  to  repeated  diffraction  of  the  waves  at  the  crack 
perimeter.  It  does  not  appear  feasible  to  carry  out  this  procedure 
analytically  since  multiple  integrals  are  involved. 

Since  it  is  clear  from  figure  (3)  that  the  condition  r>r*(t)  is 
equivalent  to  the  condition  r>r*(x),  it  follows  that  formula  (30)  de- 
fines the  stress  ahead  of  the  actual  crack  tip.  Moreover,  formula 
(30)  gives  a singularity  of  the  order  1/2  at  the  actual  crack  tip  loca- 
tion. That  this  is  so  follows  from  a Lagrange  expansion  of  the  quantity 
r-r*(r),  which  is  of  the  form  [4] 

r-r*(t)  = r-r*(t)  + I A '^*(t)  p'^(t)  , (31) 

n=l 


^E.  T.  WHITTAKER  and  G.  N.  WATSON  (1927)  Modern  Analysis,  Fourth 
Edition,  Cambridge  University  Press. 
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where 


r*(t)  = ^ , and  p(t)  = {[r-r*(t)]^  + 
ut 

On  z=o  (31 ) gives 


r-r*(x)  = !"  r + terms  of  order  . 

1 -r*(t } 

Hence  as  r^r*(t)  the  term  [r-r*(x)]"^^^  in  formula  (30)  may  be  re- 

[-1/2 


(32) 


placed  by  the  term 


r-r*(t) 
1 -r*(t) 


, so  that  in  the  neighborhood  of 


the  crack  perimeter 

k 

a(r,t)  - Ti/r[r-r*(t)J  , (33) 

where  k is  the  coefficient  of  stress  intensity 


k 


’''■*(*>  |r*(t)-t| 


f(r‘ ,t-r*(t)+r' ) 


/P~ 

/r*(t)-r‘ 


dr' 


(34) 


With  the  crack  plane  stress  completely  determined  (at  least 
for  the  time  interval  previously  noted),  it  is  now  a relatively  simple 
matter  to  obtain  the  crack  displacement.  Thus  reverting  to  the  co- 
ordinates ?=r-t,  n=r+t,  letting  w(^,n)  = v(r,t),  and  employing  the 
previous  designation  for  the  stress  component  in  formula  (17)  gives 
for  z=o 
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-2ir(^+n)  w(^,n)  = 

K 


+/ 


,a 


?*(n)  (^'-^)^/2(^+^.)1/2 


+ r T(c',n') 

n* 


N(^.'^.n'.n)  dn 


(35) 


where 


G(^' ,^,n)  = /^,  g(^' ,n' ) 


(36) 


In  (35)  £|*(ri)  denotes  the  co-ordinate  of  the  intersection  of  the 


from  (19)  the  bracketed  quantity  within  the  second  integral  vanishes. 
Hence  rewriting  (35)  in  terms  of  the  original  variables  yields 


where  region  R-|  is  shown  in  figure  (4).  Formula  (37)  represents  the 
displacement  for  general  dynamic  loading  of  a circular  crack  of  initial 
radius  a.  It  is  of  interest  to  note  that  the  solution  is  independent 
of  the  stress  ahead  of  the  crack  perimeter.  Expressions  (30)  and 
(37)  will  now  be  utilized  to  study  results  for  two  particular  cases. 


line  n'=n  and  the  crack  locus  n'=n*(E')»  i-e.  E*(n)  = . Now 


(37) 
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Figure  4.  Diagram  for  Crack  Plane  Displacement  So 


Case  1.  Initial  Value  Problem. 
In  formula  (30)  put 


f(r,t)  = 0 , o<r<a 

(0) 

= - 0{r)  , a<r<r*(t) 


(38) 


(0) 


where  a(r)  is  some  initial  state  of  stress.  The  resulting  expression 
represents  the  perturbed  stress  arising  from  the  motion  of  the  crack 
perimeter.  As  an  illustration  let 


Expression  (39)  is  the  result  given  by  Sneddon  [5]  for  the  stress  on 
the  crack  plane  z=o  under  conditions  of  uniform  stress  S at  infinity. 
Introducing  (39)  into  (30)  yields 


Since  the  initial  stress  on  the  freshly  formed  portion  of  the  crack 
surface  has  been  annulled,  then,  by  superposition,  (40)  combined  with 
(39)  gives  the  actual  stress  ahead  of  the  crack  perimeter  for  the  case 
of  uniform  stress  S at  infinity.  The  integral  may  be  expressed  in 
terms  of  known  elliptic  integrals;  these  expressions  are  given  in  the 
appendix. 

In  the  limit  as  r->-r*(t)  it  is  found  that  the  stress  intensity 

^I.  N.  SNEDDON  and  M.  LOWENGRUB  (1969)  Crack  Problems  in  the  Classical 
Theory  of  Elasticity,  Wiley  and  Sons. 


(39) 


(40) 


factor  for  this  case  is 


k 


S 

2 


r*(t) 


jr*(t) 


2r'^-a^)  dr' 

/r' (r' ^-a^)[r*(t)-r' ] 


(41) 


Explicit  expressions  for  the  integral  term  may  be  obtained  from  those 
given  in  the  appendix  by  putting  r=r*(x)  = r*(t). 

The  perturbed  crack  displacement  for  this  case  follows  by  substi- 
tuting (38)  into  (37).  Reference  to  figure  (4)  shows  that  for  o<r-t<a 
the  result  may  be  expressed  as 

yv(r,t)  = a?^)dr'  . cot  <j>  dt'  , (42) 

irr  •'a  T-r*(T)+r' 


(o) 

wherein  a(r')  is  given  by  (39)  and  t is  given  implicitly  by 

t-T  = r*(T)-r.  With  the  aid  of  tables  [3]  it  is  found  that  (42)  may 

be  written  as 

yv(r.t)  = -i-  f 'a  (r')  dr'  I(r',r,t)  , (43) 

irr  a 

where 

2 1 2 

I(r',r,t)  = F(ic,q)  - (r+r')  E{K,q) 


+ ! ([(r+r')V]  [w^-(r-r')  ])’'^  . 
w 

in  which  F(tc,q)  and  E(<,q)  are  elliptic  integrals  of  the  first  and 
second  kind,  respectively,  with 


K = sin 


r+r'  pW^-(r' -r)^-|l/2, 
w '■  4rr'  ^ 


q = 


(4rr')l/2 

r+r' 
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and  w = t-T+r*(T)-r' . For  o<t-r<a  identical  results  are  found,  hence 
(43)  is  valid  for  o<|t-rl<a. 

Case  2:  Suddenly  Appearing  Crack 

For  this  case  let  f(r,t)  = -S.  Then  formula  (30)  gives 


o(r,t)  = 


.i/r[r-r*(T)]  |r-tl 
The  result  (44),  superposed  with  a constant  stress  S,  gives  the  solu- 
tion for  a suddenly  appearing  crack  under  uniform  stress  S at  infinity. 
Here  it  is  assumed  that  the  crack  opens  instantaneously  to  radius  a 
and  thereafter  is  located  at  r=r*(t).  For  this  case  the  integral  may 
be  evaluated  explicitly;  the  result  is 


S 

o(r,t)  = /|t-r|  Lr*(T)-|t-r| J 

7r/r[r-r*(x)] 


+ •-  [2r-r*(x)]  sin"^  — — /|t-r|  [r*(x)-|t-r|] 
2 r*(x) 


- /r[r-r*(x)]  sin 


2 /r|t-r[[r*(x)-|t-r|][r-r*(x)] 


r*(x)  (r-|t-r|) 


(45) 

As  an  example,  if  the  crack  speed  is  a linear  function  of  time,  say 
r*(t)  = a+3t,  with  3<1 , then  since  t-x  = r-r*(x)  it  follows  that 


t-r+a 

^ = TT 


Putting  this  result  in  (45)  gives  for  r-r*(x)<t<r 
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/ X S ; 77 T7 T ^ (2-3)r-a-3t 

a(r,t)  = { /(r-t)(t-r+a)  + 5 — 


7r/r(r-a-3t) 


/1-3 


. -1 
sin 


V /(r-t)(t-r+a) 
a+3(t-r) 


/r(r-a-3TT 

/T3 


sin 


-1 


t[a+3(t-r)] 


/r(r-t)(t-r+a)(r-a-3t) 


(46) 


The  crack  surface  displacement  may  likewise  be  obtained  explicitly. 
Thus  substituting  into  (37)  and  referring  to  figure  (4)  gives  for 
o<r-t<a 


pv(r,t)  = - /^  dt'  cot  (j)  dr'  , r<r*(t)  . (47) 

Tfr  0 r+t'-t 

Then  using  tables  [3]  gives 


yv(r,t)  = - /^  J(r,s,t)  ds  , r<r*(t)  . (48) 

' ' irr  ■'0 

where 

J(r,s,t)  = (r-s)F(K,q)  + (r+s)E(K,q)  - - /[(r+s)^-X^][X^-r-s)^]  , 
in  which 


. -1  r+s  A^-(r-s) 

K = sin  — -7r= — — 

and  X = r*(x)  - |t-s-T| 

It  is  of  interest  to  note 


that  if  r-a<t<r+a  the  results  simplify 


^ 
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even  further,  since  for  this  time  interval  a point  on  the  crack  surface 
has  not  yet  experienced  the  effect  of  the  waves  from  the  crack  perim- 
eter. For  this  case  formula  (48)  is  replaced  by 

yv(r,t)  = — ds  cot  4>  dr'  , (49) 

irr  0 r-s 

which,  upon  evaluation  of  the  inner  integral,  gives 

vv(r,t)  = I /‘  [(r-s)K(!^)  . (r+s)E(^)]  ds.  (50) 

where  K and  E are  complete  elliptic  integrals.  Now  substituting 
s=rw  and  applying  the  transformation  formulae 

= (Hw)K(w) 

' ils 

gives 

v(r,t)  = ^ E(w)  dw  , r-a<t<r+a  . (51) 

The  final  result  follows  by  using 


(52) 

The  resulting  expression  is  equivalent  to  that  which  would  be  obtained 
on  the  surface  of  a half-space  loaded  by  a step  function  load  of 
magnitude  S. 


r^''  E(w)  dw  = 2 p 


, " [(2n)l]^(t/r)^" 

n=l  (2ntl)2''''(n!)^ 
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DISCUSSION 


The  analytical  results  contained  herein  may  be  used  as  a basis 
for  studying  the  physical  aspects  of  the  dynamic  torsional  fracturing 
process.  As  demonstrated  by  Kostrov  [2]  in  the  case  of  anti -plane 
strain,  the  crack  propagation  speed  may  be  made  definite  by  applying 
the  physical  hypothesis  that  the  work  done  in  the  rupture  process 
depends  only  on  the  crack  speed.  This  hypothesis,  together  with 
Barenblatt's  cohesion  model  leads  to  a differential  equation  governing 
the  crack  perimeter  locus. 

Another  useful  extension  of  the  results  of  this  report  is  to  incor- 
porate the  effects  of  a plastic  zone  near  the  crack  border.  This  may 
be  done  by  using  Dugdale's  model  as  shown  by  Achenbach  [6]  for  the 
case  of  anti-plane  strain.  These  considerations,  it  is  felt,  are  more 
appropriately  left  as  subjects  of  future  investigations. 


2 
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APPENDIX 

The  integral  term  in  equation  (40)  may  be  reduced  as  follows: 


Let 


„ (2x^-a^)  j 

A = J - dx  , r b , 

® A(x‘^-a‘^)  r-x 


(Al) 


where  b = r*(x).  Changing  the  integration  variable  to  u = r-x  leads  to 

(A2) 


A = -a"I  + 2J 


where 


I = I1-3I2 


J = 


(A3) 


(A4) 


du 


I,  = r =zzz 

3 /(r-u)(a-u)(u-3)(c-u) 


(A5) 


du 


1 = r— 

2 3 u/(r-u)(  -u)(u-  )(c-u) 


(A6) 


a *^r-u^  du 


3 /(a-u)(c-u)(u-3) 


(A7) 


Jo  = /; 


a 


/r-u^  du. 


^ 3 u/(a-u) (c-u) (u-3) 

In  (A3)  - (A8),  3=r-a,  a=r-b,  c=r+a,  so  that  c>r>a>3  . Now  from 
tables  [3]  it  is  found  that 


(A8) 


3 
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Ac-yjw-&) 

I,  = /(c-a)(r-B)  K( ] 

' ^lc-a){r-B) 


(A9) 


where  K is  the  complete  elliptic  integral  of  the  1st  kind.  Hence 
replacing  the  original  variables  gives 

2 


^1  = 


/r*(x)-a 

K(-=:r)  • 


/2ar*(x)  »^2r*U) 


(AlO) 


I = 2[a-»-r*(x)] 

^ (r+a)[r-r*(x)]/2ar*(x)  2a[r-r*(x)]  2r*(x) 

(All) 

when  II  is  the  elliptic  integral  of  the  third  kind.  To  evaluate  J, 

u “ 6 

make  the  further  change  of  integration  variable  z = , so  that 

(A7)  becomes 


j = /(r-B)^  r^  /(1-6z)3  dz 
^ /c^  "^0  /z(l-z)(l-rz) 


(Al2) 


in  which  6 = , and  y = . Formula  (A12),  as  given  in 

tables  [3],  may  be  expressed  as 


, . /[r*(x)]^  . .1  31, 

1 ■ /fPFTET  ’''^1(2’-  2’  2-  '■ 

where  F-j  is  a hypergeometric  function, 
evaluated  by  writing 

J2  = r - J22  . 


r*(x)-a  r*(x)-a>> 

r*(x)  ’ r*(x)+a-' 

Finally,  integral  J2 


is 


(A13) 


(A14) 


3 
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In  equation  (A14) 

in  which  I2  and  I-|  are  given  by  (AlO)  and  (All),  and 


J 


22 


/r^  du 
/(a-u)(c-u)(u-B) 


(A15) 


(A16) 


The  integral  J22  is  given  directly  by  tables  [3];  the  result  in  terms 
of  the  original  variables  is 


1 = n(ir/?  /r*'(T)-a-> 

22  /r*(t)  * ”r*T^T~  ’ ,^zr*{x) 

Substituting  (A9)-(A17)  into  (A2)-(A4)  gives  the  result  for  (Al). 


(A17) 


3 
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